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Abstract 

An optimization-based methodology for linear control system design is outlined by con- 
sidering the design of a controller for a UH-60 rotorcraft in hover. A wide range of design 
specifications is accounted for: internal stability, decoupling between longitudinal and lat- 
eral motions, handling qualities and rejection of wind gusts, while taking into account physi- 
cal limitations in the swashplate displacements and rates of displacement. The methodology 
crucially relies on user-machine interaction for tradeoff exploration. 

1. Introduction 

Airframe concepts for future rotorcraft such as the light attack helicopter (LHX) include 
high effective hinge offset rotors and multiple control effectors to maximize maneuverability. 
This and other possible features lead to a high level of bare airframe instability, control cross 
coupling, and control redundancy. At the same time there is a need to include explicitly in 
a design methodology numerous specifications on control limits, actuator capabilities, cross 
coupling limits, and handling qualities. The classical control techniques which have been the 
primary tools of the industry are not well suited for such complex combinations. There is a 
great need for improved techniques. Numerous modern control methodologies have recently 
been proposed for application to these types of problems(e.g., [1-3]). However these studies 
have not adequately accounted for the many practical implementation problems of rotorcraft 
[4]. Key concerns in evaluating prospective modern control designs are the interaction of the 
rotor system dynamics [5] and the explicit inclusion of the complex design specifications( Mil- 
Specs [6]) in the design process. 

A methodology is proposed that can account for various types of concurrent specifica- 
tions: stability, decoupling between longitudinal and lateral motions, handling qualities, 
and physical limitations of the swashplate motions. This is achieved by synergistic use of 
analytical techniques (Q-parametrization of all stabilizing controllers, transfer function in- 
terpolation) and advanced numerical optimization techniques. The proposed methodology 
includes a two-parameter controller with separate consideration of the feedforward and the 
feedback portions. Preliminary results on the input/output part of the design were reported 
in [7], where a 21-state model was used for the aircraft. Here overall results are reported on 
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a 37-state model generated by the UM-GenHel [8] nonlinear simulator. Nonlinear validation 
is also presented. 

A central tool in this study was the interactive optimization-based package Console 
[9,10]. CONSOLE handles multiple (competing) objective functions, and constraints can be 
“soft” (moderate violations are allowed) or “hard”. Specifications can be “functional”, i.e., 
involve an entire time- or frequency-response rather than a single time or frequency point. 
For each specification, the designer provides a “good value” and a “bad value” (a “good 
curve” and a “bad curve” in case of a functional specification) and these can be interactively 
adjusted as tradeoffs are identified. 

In Section 2 below, we briefly describe a simplified model of a rotorcraft in hover (based 
on the UH-60). In Section 3, a set of design specifications is proposed. In Section 4, 
the design methodology is outlined. In Section 5, some results are presented, including 
validation on the nonlinear model. Section 6 is devoted to some final remarks. 

2. A Linearized Model for an UH-60 in Hover 

The model we used for controller design, denoted by Pq(s), is a 39-state linear model 
generated from UM-GenHel. UM-GenHel, developed at the Aerospace Engineering Depart- 
ment of the University of Maryland, is a nonlinear flight simulation program for helicopters 
[8]. It models the high order dynamics (such as main rotor blade motions and main rotor 
inflow) of the helicopter, and can generate linearized models of the helicopter at various 
flight conditions. The last six states of this linearized model correspond to hidden modes 
associated with the engine dynamics. These six states were removed from the model, leav- 
ing a 33-state minimal realization of Pq($). This realization has one pole at zero and a pair 
of real poles in the open right-half plane. It is non-minimum phase and strictly proper. 
In each of the four channels, the control system actuator is modeled by a first order Pade 
approximation corresponding to a 0.050 second delay. This yields a final count of 37 states. 

An overall block diagram of the closed-loop system is represented in Figure 1. The 
four variables of the input 6 S = (6 s< t>, <5 5 <9, S 5C , 6 S ^) of Pq(s) represent respectively the lateral, 
longitudinal, and collective displacements of the swashplate, and the position of the tail 
rotor actuator. The output is y = (u y v y w,p,q,r, 4>,0, 0), the variables of which stand for 
longitudinal velocity (u), lateral velocity (u), vertical velocity (te), roll rate (p), pitch rate 
( q ), yaw rate (r), roll angle (0), pitch angle (0), and yaw angle (0). The delay block D e (s) 
has 4 states, 4 inputs, and 4 outputs. The command input S — (<5^, <5$, 6 C , 6^) consists of 



Figure 1: Control system configuration 


380 






roll command (6^), pitch command (£*), collective command (<5 C ), and yaw command (6$). 
The rotorcraft is to be controlled by a two-parameter controller C(s). 

3, Design Specifications 

A wide range of specifications, in both the time- and frequency domains, are to be 
satisfied. First, the closed-loop system must be internally stable. Second, to the extent 
possible, it is desired that the various longitudinal and lateral modes be decoupled and 
approximate specified step responses (handling quality specifications). Specifically, the 
pitch command should mostly affect the pitch angle (and pitch rate) and the longitudinal 
velocity; the collective command should mostly affect vertical velocity; the roll command 
should mostly affect the roll angle (and roll rate) and the lateral velocity; and the yaw 
command should affect mostly the yaw angle (and yaw rate). The “diagonal” responses 
should exhibit desirable characteristics as given in [6]. Third, the displacement and rate 
of displacement of the swashplate may not violate some physical limitations. Lastly, the 
closed-loop system should exhibit good performance in gust rejection. The mathematical 
formulation of these specifications is described below. (Note: all inputs are expressed in 
inches of stick.) 

Spec 1 The closed-loop system is internally stable. 

Spec 2 When the input 6 = [0 5 0 0] T : 

Spec 2.1 Pitch response 6(t) should be between the two curves C\(i) and C 2(0 of Figure 2. 



steady state (s.s.) = 0.5 rad. 
overshoot = 35% settling = 5% 


Figure 2: Boundaries for pitch response 

Spec 2.2 Decoupling: (i) <f>(t), rft(t) should have the absolute values less than 5% of 0(5) for 
t £ [0,5]; (ii) v(t), w(t) should have the absolute values less than 5% of u(5) for t £ [0,5]. 
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Figure 3: Handling quality specification 

Spec 2.3 Physical limitations on swashplate displacements: the absolute values of S 3 g(t) + 

MO. MO + MO. MO + MO. MO. MO. MO. and MO should al] be less than 10 

inches for / G [0,5]. 

Spec 2.4 Limitations on the velocities of the swashplate motion: the absolute values of 
M0 + M0>M0 + M0,M0 + M0,M0A*(0>M0>and MO should all be less 
than 10 inches per second for t G [0,5]. 

Spec 3 Let Hg be the transfer function from 6g to 6. Suppose ubw (bandwidth) and t v 
(phase delay) are defined based on the Bode plot (phase part) of Hg , such that ubw is the 
frequency corresponding to 45 degree phase margin, and phase delay t p is defined as 

_ $2u>ig 0 + 180 

Tp ~ X 2u xm 

where uqso is the frequency at which the phase angle attains —180 degree, and $2 wi 80 ’ s 
the phase angle at frequency 2u;j8o (see [4]). Then the graph of (mbw,t p ) should lie within 
region Level 1 in Figure 3. 

Spec 4 When the input 6 = [5 0 0 0] T : 

Spec 4.1 Roll response <f>(t) should lie between the two curves Ci(t) and C 2 (i) of Figure 2, 
with s.s.= 1 rad. 

Spec 4.2 Decoupling: (i) 0(t), ip(t) should have the absolute values less than 5% of <j>( 5) for 
t G [0,5]; (ii) u(t), w(t) should have the absolute values less than 5% of v(5) for t G [0,5]. 
Spec 4.3 Physical limitations on swashplate displacements: same as spec 2.3. 

Spec 4.4 Limitations on the speed of swashplate motion: same as spec 2.4. 

Spec 5 Same as Spec 3 with Hg replaced by H the transfer function from 6$ to </>. 
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Spec 6 When input = [0 0 0 2.5] : _ 

Spec 6.1 Yaw rate r(t) should lie between the two curves Ci(t) and C2(0 of Figure 2, with 

s.s.= 0.4 rad/sec. 

Spec 6.2 Decoupling: p(t ), q(t) should have the absolute values less than 5/o of r(5) lor 
t € [0, 5]. 

Spec 6.3 Physical limitations on swashplate displacements: same as spec 2.3. 

Spec 6.4 Limitations on the speed of swashplate motion: same as spec 2.4. 


Spec 7 Same as Spec 3 with H e replaced by H+, the transfer function from to 0. 

Spec 8 When input = [ 0 0 5 0 ] : 

Spec 8.1 Vertical velocity w(t) should satisfy 


' bKe~ TS ' 

(f \ / < 1 1 

■ 5Ke~ Ta ' 

_s(Ts+ 1)_ 

^ Z J ^ W \ t J _ * f 

s(Ts + 1) 


(0 + 2 


for all t € [0,5], where K = 10, r = 0.2, T - 5, and C~ l denotes the inverse Laplace 
transform. 

g pec g,2 Decoupling: u(«), v(t) should have the absolute values less than 5% of u>(5) f° r 
f G [0,5]. 

Spec 8.3 Physical limitations on swashplate displacements: same as spec 2.3. 

Spec 8.4 Limitations on the speed of swashplate motion: same as' spec 2.4. 

Spec 9 A wind gust, which is modeled by the step response of 


g(s ) = 7~7T^ (6 ust model) (1) 

(s + ir 

is to be injected at the output nodes <f > , 8, and 0 in turn. 

Spec 9.1 The closed-loop transient response at the corresponding node, when the wind gust 
is injected, should be kept within the envelope of the step response of 


0.155s 

haiiowfs) - 5 2 + 1.77 S + 0.462' 


( 2 ) 


Step responses of g(s) and h a \] ovl (s) are shown in Figure 4. 

Spec 9.2 Physical limitations on swashplate displacements, same as spec 2.3, but time in- 
terval is changed to [0,oo). 

Spec 9.3 Limitations on the speed of swashplate motion, same as spec 2.4, but time interval 
is changed to [0, oo). 


We remark that, given the highly coupled dynamics of the plant (P 0 {s)), these speci- 
fications are very stringent, and to our knowledge no controller has been built so far that 
satisfies all these specifications. 

The asymmetry of the rotorcraft places fundamental limitations on the amount of de- 
coupling that can be achieved between the various channels. The most striking instance is 
the impossibility of maintaining a constant nonzero pitch angle without steady state bank- 
ing. For the linear model used in this study, it can be verified (see [11,12] for details) that 
the minimum achievable ratio between steady state roll and pitch angles is slightlj o\er 
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Figure 4: Gust model and allowable closed-loop response 


20%. Note however that this does not preclude the possibility of meeting specification 2.2 
as it requires significant decoupling over the first 5 seconds only. Yet, this limitation ought 
to be kept in mind when interactively exploring possible tradeoff designs. 

4. Methodology 

Design of a two-parameter controller proceeds as follows [13]. First the overall aircraft 
transfer matrix (rotor + airframe dynamics + delay) P(s) = P 0 (s)D e (s) is factorized over 
the ring of stable transfer functions, i.e., 

P(s) = N(s)D-\s ) = D~ 1 (s)N(s ) 

where N($)> D(s), N(s) and D(s) are stable transfer matrices, (N(s), D(s)) are right 
coprime and (jV(s), D(s)) are left coprime. Next let A"(s) and Y(s) be stable transfer 
matrices satisfying the Bezout identity 

X(s)N(s) + Y(s)D(s) = / 

Then all stabilizing controllers can be obtained according to the block diagram of Figure 5 
where both Q(s) and R(s ) range over the class of all stable transfer matrices of appropriate 
dimension (with the condition that det(Y(5) - R(s)N(s)) does not vanish identically). 



Figure 5: Two-parameter controller 
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It is easily verified that the overall transfer function T(s) from S to y is given by 

T{s) = N(s)Q(s) (4.1) 


while the transfer function from 6 to the swashplate displacements 6 S by 

S(s) = D e (s)D(s)Q(s) (4.2) 

Thus both are independent of the R parameter. On the other hand the transfer function 
from disturbance inputs d e and d y to those same outputs y and 6 a is given by 

r Z>(*)(y(«) - R(s)N(s)) -D(s)(X(s) + R(s)D(s)) ■ 

|7V(s)(r(s) - R(s)N(s)) -N(s)(X(s) + R{s)D(s)) + /. 

which is independent of the Q parameter. By adopting the two-parameter controller struc- 
ture, the control designer can then separate the design problem into two independent sub- 
problems of input/output performance (to be solved by choosing a ’’good” Q: Q-design) 
and of disturbance rejection (to be solved by choosing a "good” R: 72-design), respectively. 
At the same time, internal stability will be automatically guaranteed if the chosen Q(s) and 
72(.s) are stable. 

Q-design 

We just saw that any stable Q(s) yields a stable T(s) and any stable T(s) corresponds 
to a stable Q(s), so that this “Q-parameterization” automatically takes care of specification 
1. Concerning the other input/output specifications, the structure of the problem allows 
another major simplification. Note that each one of the specifications 2-8 involved only one 
of the four input channels and (see (4.1-2)) each column of Q(s ) affects the I/O 
transfer function corresponding to exactly one input channel (i.e., exactly one column of 
Q(s)). As a result the problem can be decomposed into four optimization problems, each 
one involving a single column of Q(s) and the specifications involving inputs corresponding 
to that column. Thus, interactive optimization will be used, for each input channel, to try 
to satisfy the corresponding input/output specifications; this will involve a single column 
of Q(s), which we will denote by q(s). 

A possible parametrization for q(s), once its McMillan degree n has been chosen, is by 
writing 

q{s) = C(sl- A)~ l B+ D (4.3) 

with _ _ _ 


'01 


V 

Pi P2 


0 

0 1 


l 

P3 Pa 

, B = 

0 

0 1 


1 

P2j-1 P2j . 


_ 0 _ 
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if n = 2 j is even, and 

0 1 1 r 1 

Pi P2 0 

•4= 0 j ■ « = 0 

P2j-l P2 } '■ 

P2j+1 J 1 - 

if n = 2j + 1 is odd, and using the pf s as well as all entries of C and D as free parameters. 
Indeed, it is clear that by constraining the p,’s to be negative, only stable Q(s) will be 
generated; it is also easy to show that the family of matrices Q(s) thus generated is open 
and dense in the set of stable matrices of McMillan degree n. 

One difficulty with the parametrization just proposed is the ensuing large number of 
design parameters, resulting in (i) large computation times and (ii) likeliness that a local 
rather than global optimum be reached. The following heuristics drastically reduces the 
number of design parameters by eliminating a large portion of “clearly nonoptimal” q(s) 
and likely retaining enough of the “good” ones to allow a “close to optimal” q(s ) to be 
identified. For given A and B (i.e., given p;’s) and given “sampling frequencies” 

(the p^s and the u>j’s will be the design parameters) corresponding C and D are selected 
by solving, in a least square sense, 

N(ju,)(C{ju,I - A)~ l B + D) = T(ju t ), i = 1 , . . . , * (4.5) 

where T(s ) is an approximate “desired” transfer function constructed from the specifica- 
tions: good “diagonal” responses and “off diagonal” responses set to zero. (The method that 
was actually used is somewhat more sophisticated in that it involves solving two successive 
linear least squares problems: see [11,12].) 

To summarize, the numerical optimization process will proceed as follows. Given values 
of the p^s and u>i’s (design parameters), q(s ) will be obtained by solving for C and D the 
linear least square problem; time and frequency responses of the closed-loop system will 
be computed (by invoking MATLAB); the values of the specifications and their gradients 
(with respect to the design parameters) will be determined; the optimization algorithm (the 
heart of Console) will use this information to select new values of the design parameters; 
and the cycle will start again until the designer decides to interrupt it. The designer can 
then relax or tighten selected specifications in order to explore various tradeoffs. 

In addition to the specifications outlined in Section 3, it may be desired to achieve or 
approach certain steady state values when unit steps are fed into individual channels. In 
view of the Final Value theorem, these values are determined by <2(0). The best value Q, 
for (J(0) can be determined, again, by solving a linear least squares problem. Let q « be a 
column of Q ,. Then, in view of (4.3), we may require that 

C(-A)-'B + D = q., 

i.e., 

D = g. + CA~ l B. 
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Substituting this expression for D in (4.5) we obtain a linear least squares problem in C only. 
Note that, in the design results reported below, there was no such steady state requirement. 

/{-design 

A major difference with the Q-design is that the columns of /{(.s) cannot be designed 
independently, as each of these columns may affect all disturbance rejection specifications. 
Thus a parametrization of the form (4.3-4), where C and D are obtained by solving a linear 
system of the form (4.6) in the least squares sense, is used for the entire matrix R(s). 

5. Results and Validation 

Based on our Q-design results, it is found that in the lateral, longitudinal, and tail 
collective channels, the swashplate rate specifications (Specs 2.4, 4.4, 6.4) are competing 
with the handling quality (based on frequency response) specifications (Specs 3, 5, 7). 
Specifically, for the handling quality indicator (w B w,r p ) to lie within region Level 1 in 
Figure 3, the swashplate rates will attain magnitude above 10 inch/sec within the first 0.1 
second. For example, Figure 6 and Figure 7 show the results of a design of the second 
column of Q (the longitudinal channel) on Specs 2.4 and 3. For this design, Spec 2.4 is 
satisfied but Spec 3 is not. Figure 8 and Figure 9 show the results of another design of the 
same column of Q on Specs 2.4 and 3. For this design, Spec 3 is satisfied but Spec 2.4 is 
not. 

Our /{-design results are that Specs 9.2 and 9.3 are satisfied, but Spec 9.1 is not. Figure 
10 shows the response of the roll angle <f>(t ) (or, respectively, the pitch angle 0(t), or the yaw 
angle when a wind-gust, modeled by the step response of 0.44s/(s + l) 2 , is injected at 

the roll angle node (or, respectively, the pitch angle node or the yaw angle node). From the 
figure, one sees that the graphs of the response do not lie within the allowable envelope given 
in Spec 9.1. Hence Spec 9.1 is not satisfied. The figure also suggests that the helicopter 
cannot act fast enough to respond to the gust, and can counteract the gust only after 0.3 
second. 

We finally choose a design of the controller for which the specifications on swashplate 
rates (Specs 2.4, 4.4, and 6.4) are satisfied but the handling quality specifications (Specs 
3, 5, and 7) are not. This controller has 152 states. It is then reduced to 82 states via the 
balanced realization method. It has been examined, by checking against the specifications, 
that the performance of the order-reduced controller on the linearized model is as good 
as the unreduced, 152-state controller. This order-reduced controller is tested on both the 
UM-GENHEL linear and nonlinear simulation program. Based on the comparison between 
the open-loop time responses (to a step input) of the nonlinear and the linearized model, it 
is found that the nonlinearity of the helicopter model is very high. 

Figure 11 shows some of the relevant time responses of the closed-loop (linear and 
nonlinear) systems, when a 5 unit (in terms of inches of stick command) step command 
input is put on the lateral channel (£,*,) of the controller. One notices from this figure that, 
after 2.5 seconds, the closed-loop responses of the nonlinear model do not match those of 
the linear model very well. However, given the high nonlinearity of the dynamics of the 
nonlinear model, this phenomenon should be expected. For example, Figure 11 shows that 
at 2.5 seconds after the control command has been issued, the roll angle of the nonlinear 
model is about 1 radian. This flight condition is very different from hover, and hence 
the dynamics of the nonlinear model at this point may be very different from that of the 
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linearized model for hover. 

One would expect the discrepancies between the response of the closed-loop nonlinear 
system and that of the linear system should be less, if the magnitude of the control input 
is lowered. Figure 12 gives the comparison between the relevant responses of the nonlinear 
and linearized model, when the magnitude of the command input is lowered to 1. This 
clearly confirms that the response of the nonlinear model matches better with that of the 
linear model than before. 

Finally, closed-loop responses of both the nonlinear and linear model, due to gust dis- 
turbance, are examined. For example, Figure 13 shows the roll angle response when gust 
is injected at the roll angle node. Notice that there are some discrepacies between the roll 
angle responses of the linear and the nonlinear model. 

6. Discussion 

The controller discussed in Section 5 is but one of the many sub-optimal controllers 
obtained when running Console. Indeed, a key advantage of an interactive package such 
as CONSOLE is that it allows the designer to explore alternative solutions by fine tuning the 
various target responses (this is especially so with Console’s not yet released graphical 
interface). 

It turns out that even the stability specification is amenable to tuning. Indeed, as 
discussed in [13], the Q-parametrization approach can be extended to generalized stability, 
i.e., confinement of closed loop poles in a more general region II of the complex plane. This 
can be achieved by (i) factoring P(s) in the ring of IT-stable (rather than Hurwitz stable) 
transfer functions and (ii) take for <2(s) and R(s) any matrices with all their poles in II. 
The latter can be accomplished by suitably modifying parameterization (4.4). E.g., if IT is 
as in Figure 14, it suffices to replace in (4.4) 2x2 blocks of the form 

0 1 
Pi P2 

0 1 ' 

(1 - y)(b - x) 2 - x 2 (yk 2 + 1) 2z. 

and let x and y vary over [(a + 6)/2, 6] and [0, 1], respectively (see [11,12] for details). 

Finally, the merits of the approach discussed in this paper should be compared to those 
of the convex optimization approach proposed by S.P. Boyd and C.H. Barratt [14], The 
main advantage of the latter is that it always yields the globally optimal design for the 
given specifications (for Q(s ) and R(s) ranging over a finite dimensional subspace of the 
space of stable transfer matrices). A crucial requirement, however, is that all specifications 
be convex (as functions of the closed loop transfer matrix). It turns out that some of the 
specifications considered in the present study (in particular handling quality specifications 
such as Spec 3) do not satisfy this requirement (see [11,12] for details). Compared to the 
Ritz parameterization used in [14], while parameterization (4.3-4) would destroy any existing 
convexity, it has the advantage to cover all (in fact, “almost all”) stabilizing controllers of 
a given degree. 
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Figure 6. Pitch command: swashplate rates 
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Figure 7. Bandwidth and phase delay for pitch channel 


Spec 2.4a, 2.4b, 2.4c: Swashplate rate 
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Figure 8. Pitch command: 
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Figure 10. Time- responses to wind gust 
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5 inches input at lateral channel 
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Figure 11. Nonlinear validation: large signals 
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Figure 12. 
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